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Heavy quark jet quenching in nuclear collisions at LHC is predicted and compared using the clas- 
sical gravity AdS/CFT correspondence and Standard Model perturbative QCD. The momentum 
independence and inverse quark mass dependence of the drag coefficient in AdS/CFT differs sub- 
stantially from the characteristic log (pr /M) /pr variation of the drag in QCD. We propose that the 
measurement of the momentum dependence of the double ratio of the nuclear modification factors 
of charm and bottom jets, Raa(Pt) / R\a{Pt) , is a robust observable that can be used to search for 
strong coupling deviations from perturbative QCD predictions. 
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Introduction. Recent discoveries at RHIC pjj EJ [3] 
have led to provocative suggestions j4] that the proper- 
ties of strongly coupled quark gluon plasmas (sQGP) pro- 
duced in ultra-relativistic nuclear collisions may be bet- 
ter approximated by string theoretic inspired AdS/CFT 
gravity-gauge theory correspondences [5] than conven- 
tional Standard Model perturbative QCD (pQCD). Four 
main classes of observables have attacted the most atten- 
tion: (1) Entropy production as probed by multiplicity 
distributions [3] , (2) Perfect Fluidity [3] as probed by col- 
lective elliptic flow observables HJ [3] , (3) Jet Quench- 
ing as probed by high-p^ hadrons jTUJ E] and leptons 
[12), H"3] . and (4) Dijet-Bulk Correlations as probed by 
two and three particle correlations [TU [TS] . 

Qualitative successes of recent AdS/CFT applications 
[3] to nuclear collision phenomenology include the ana- 
lytic account for (1) the surprisingly small (~ 3/4) drop 
[16] of the entropy density in lattice QCD calculations 
relative to Stefan-Boltzmann, (2) the order of magnitude 
reduction of the viscosity to entropy ratio t]/s predicted 
relative to pQCD needed to explain the seemingly near 
perfect fluid flow of the sQGP observed at RHIC, (3) the 
surprising large stopping power of high transverse mo- 
menta heavy quarks as inferred from the quenching and 
elliptic flow of leptonic decay fragments of heavy quark 
jets, and (4) the possible occurrence of conical "Mach" 
wave-like correlations of hadrons associated with jets. 

While quantitative and systematic comparisons of 
AdS/CFT gravity dual models with nuclear collision 
data are still incomplete and while the conjectured dou- 
ble Type IIB string theory<->conformal Supersymmetric 
Yang-Mills (SYM) gauge theory<-s-non-conformal, non- 
supersymmetric QCD correspondence remains under de- 
bate (see, e.g. [T7]), the current successes provide strong 
motivation to seek more sensitive experimental tests that 
could help guide the development of such novel theoreti- 
cal approaches. 

The aim of this Letter is to propose a robust new test 
that could reveal possible strong coupling deviations from 



pQCD predictions as suggested by a specific version of 
the AdS/CFT correspondence. This test involves the nu- 
clear modification of identified heavy quark jets produced 
in central Pb+Pb reactions at 5.5 ATeV at LHC. Similar 
tests can be performed at RHIC, after future detector up- 
grades, and will be reported elsewhere. Specifically, we 
propose that the double ratio of identified charm and bot- 
tom jet nuclear modification factors R^ a (pt) can easily 
distinguish between a wide class of pQCD models and a 
class of gravity dual models [T5J EH 1201 ES [22 of heavy 
quark dynamics. The current failure of pQCD models 
to account quantitatively for the recent RHIC data from 
STAR [12] and PHENIX [13^ on the nonphotonic electron 
spectrum provides additional strong motivation to focus 
on heavy quark jet observables. Unlike for light quark 
and gluon jet observables, where pQCD predictions were 
found to be remarkably quantitative |llj . heavy quark 
jet quenching, especially as inferred indirectly for bot- 
tom quarks, appears to be significantly underpredicted 
[533 HI H5]. However, the pQCD predictions can not 
yet be falsified at RHIC because of (1) remaining un- 
certainty in the nuclear production ratio of bottom to 
charm quarks and (2) current controversy over the rel- 
ative magnitude of elastic versus radiative loss channels 
[213 12S]- Recent attempts to reconcile pQCD predictions 
with the data by adding nonperturbative hadronic final 
state interaction effects can be found in J57J [2B] . 

The AdS/CFT correspondence has so far been applied 
to QCD jet physics in three different ways. The first 
calculates the Wilson line correlator needed to predict 
the partonic radiative transport coefficient q [29l [30] . 
The second predicts the heavy quark diffusion coeffi- 
cient D |31j used as an input parameter in a relativistic 
Langevin model with drag deduced from the Einstein re- 
lation [32] ■ The third computes directly the heavy quark 
drag coefficient via a specific classical string configuration 
p~9| l20l l2Tj . We note that the first two approaches con- 
form more closely to the original spirit of the AdS/CFT 
conjecture relating quantum SYM correlation functions 
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to the asymptotic behavior of the classical supergrav- 
ity (SUGRA) correlations. Unfortunately these corre- 
spondences are hard to interpret in terms of gauge the- 
ory energy loss mechanisms since infinitely coupled SYM 
does not support a familiar quasiparticlc basis similar to 
gluon and quark degrees of freedom in QCD. In contrast, 
the third model, while more easily interpretable, requires 
a stronger form of the AdS/CFT correspondence. All 
three approaches remain under active debate (see, e.g. 

unsniiii)- 

We focus in this Letter on the third proposed AdS/ 
CFT application that involves the most direct string the- 
oretic inspired gravity "realization" of heavy quark dy- 
namics [HJ [HI [20] . A heavy quark in the fundamental 
representation is a bent Nambu-Goto string with one end 
attached to a probe brane and that trails back above the 
horizon of a D3 black brane representing the uniform 
strongly coupled SYM plasma heat bath. This geometry 
maps the drag force problem into a modern string the- 
oretic version of the old 1696 Brachistochrone problem 
that yields a remarkable, simple analytic solution for the 
string shape and momentum loss per unit time. 

AdS/CFT compared to pQCD Exploiting this AdS/ 
CFT correspondence, the drag coefficient for a massive 
quark moving through a strongly-coupled SYM plasma 
in the A = g 2 SY M N c > 1, iV c > 1, M Q > T* limit is 
given in [H [Ml IS] as 



dpT 
dt 
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where T* is the temperature of the SYM plasma as fixed 
by the Hawking temperature of the dual D3 black brane. 
Issues related to the relaxation of the strong assumptions 
made in deriving, and the momentum limitation of the 
applicability of, Eq. will be discussed later in the text. 

Applying Eq. to LHC requires an additional pro- 
posal that maps QCD temperatures and couplings to 
the SYM world and its SUGRA dual. The "obvious" 
first prescription |35j is to take gsYM — 9s constant, 
T* = T QCD , and N c = 3. However it was suggested in 
[35] that a more physical "alternative" might be to equate 
energy densities, giving T* = T < 3 CI) /3 1 / 4 , and to fit the 
coupling A = g 2 SY M^c ~ 5.5 in order to reproduce the 
static quark-antiquark forces calculated via lattice QCD. 

The string theoretic result for the diffusion coefficient 
used in the Langevin model is D/2ttT* = [S]. 
This illustrates well the problem of connecting the T* 
and A of SYM to "our" QCD world. Using the "obvious" 
prescription with a s — .3, N c = 3, one finds D/2ttT ~ 
1.2. However, D/2irT = 3 was claimed in [T3J [3T] to fit 
PHENIX data somewhat better. Note that D/2ttT = 3 
requires an unnaturally small a s ~ 0.05 that is very far 
from the assumed A> 1 't Hooft limit. 

We proceed by computing the nuclear modification fac- 
tors, neglecting initial state shadowing or saturation ef- 
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FIG. 1: (Color Online) Raa{Pt) and Raa{pt) predicted for 
central Pb+Pb at LHC comparing AdS/CFT Eq. and 
pQCD using the WHDG model [25] convolving elastic and in- 
elastic parton energy loss. Possible initial gluon rapidity den- 
sities at LHC are given by dNJdy = 1750, from a PHOBOS 
PI55] extrapolation, or dN g /dy = 2900, from the KLN model 
of the color glass condensate (CGC) [39 . The top two curves 
from pQCD increase with pr while the bottom two curves 
from AdS/CFT slowly decrease with p T . The AdS/CFT pa- 
rameters here were found using the "obvious" prescription 
with otsYM = -05, t = 1 fm/c, giving D/2-kT — 3 (abbrevi- 
ated to D — 3 in the figure). Similar trends were seen for the 
other input parameter possibilities discussed in the text. 



fects. In order to correctly deconvolute such effects from 
the final state effects that we compute below, it will be 
necessary to measure nuclear modification factors in p+A 
as a function of (u,Pt) at LHC just as d+A was the crit- 
ical control experiment [T] at RHIC [2]. 

Final state suppression of high-py jets due to a frac- 
tional energy loss e, p^ T = (1 — e)p l T , can be com- 
puted knowing the Q-flavor dependent spectral indices 

"' lo § f^fe) from PQCD or directly 
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from p + p — ► Q + X data. The nuclear modification fac- 
tor is then R® a (pt) =< (1 — e)" Q >, where the average 
is over the distribution P(e; Mq,px, I) that depends in 
general on the quark mass, pt, and the path length £ of 
the jet through the sQGP. As in [55] we average over jets 
produced according to the binary distribution geometry 
and compute I through a participant transverse density 
distribution taking into account the nuclear diffuseness. 
Given dN g /dy of produced gluons, the temperature is 
computed assuming isentropic Bjorken ID Hubble flow. 
As emphasized in [25 , detailed geometric path length 
averaging plays a crucial role in allowing consistency be- 
tween 7r°,?7 and heavy quark quenching in pQCD. 

For AdS/CFT drag, Eq. gives the average frac- 
tional energy loss as e = 1 — exp(— /j,q£). Energy loss is 
assumed to start at thermalization, tq ~ 0.6 — 1.0 fm/c, 
and stops when the confinement temperature, T c ~ 160 
MeV, is reached. The exponentiated T 2 dependence in 
fiQ leads to a significant sensitivity to the opacity of the 
medium, as well as to To and T c . 

To understand the generic qualitative features of our 
numerical results it is instructive to consider the simplest 
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case of a geometric path average over a static, finite, uni- 
form plasma of thickness L; then 



R Q aa(pt) = 



1 



(2) 



where the px dependence is carried entirely by the spec- 
tral index uq{pt). Raa can be interpreted for L^> £q = 
l/(nQHo) as the fraction £q/L of the Q jets that escape 
unstopped from the strongly coupled plasma within the 
AdS/CFT approximation. 
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FIG. 2: The double ratio of R%_a(pt) to R b AA (pr) predictions 
for LHC using Eq. Q for AdS/CFT and WHDG [25] for 
pQCD with a wide range of input parameters. The generic 
difference between the pQCD results tending to unity con- 
trasted to the much smaller and nearly pr-independent results 
from AdS/CFT can be easily distinguished at LHC. 

Two implementations of pQCD energy loss are used in 
this paper. The first is the full WHDG model convolving 
fluctuating elastic and inelastic loss with fluctuating path 
geometry [25]. The second restricts WHDG to include 
only radiative loss in order to facilitate comparison to 
[30] . Note that when realistic nuclear geometries with 
Bjorken expansion are used, the "fragility" of Raa for 
large q reported in |36] is absent in both implementations 
of WHDG. 

Unlike the AdS/CFT dynamics, pQCD predicts 
[2"31 [Ml 125] that the average energy loss fraction 
in a static uniform plasma is approximately e sa 
nL 2 q\og{pT /Mq) /pt, with n a proportionality constant 
and q = fijj/Xg. The most important feature in pQCD 
relative to AdS/CFT is that c p qcd — * asymptotically 
at high-pr while eAds remains constant. «-q(pt) is a 
slowly increasing function of momentum; thus R V aa CD 
increases with pr whereas Raa decreases. This generic 
difference can be observed in Fig. [I] which shows repre- 
sentative predictions from the full numerical calculations 
of charm and bottom Raa(pt) a t LHC. 

Double Ratio of charm to bottom Raa disadvantage 
of the Rj[ a (pt) observable alone is that its normaliza- 
tion and slow pt dependence can be fit with different 
model assumptions compensated by using very different 



medium parameters. In particular, high value extrapola- 
tions of the q parameter proposed in [26] could simulate 
the flat pr independent prediction from AdS/CFT. 

We propose to use the double ratio of charm to bot- 
tom Raa to amplify the observable difference between 
the mass and pr dependencies of the AdS/CFT drag 
and pQCD-inspired energy loss models. One can see in 
Fig. [2] that not only are most overall normalization dif- 
ferences canceled, but also that the curves remarkably 
bunch to either AdS/CFT-like or pQCD-like generic re- 
sults regardless of the input parameters used. 

The numerical value of R cb shown in Fig. [2] for 
AdS/CFT can be roughly understood analytically from 
Eq. Q as, 



R cb 



AdS 



M c n\ PT ) 
M b n c (p T ) 



Mc 



0.26, 



(3) 



where in this approximation all A, T* , L, and n c (pr) « 
n b(PT) dependences drop out. 

The pQCD trend in Fig. [2] can be understood qualita- 
tively from the expected behavior of £ p qcd noted above 
giving (with n c w nj, = n) 



n-pQCD 



Pcb 
PT ' 



(4) 



where p cb — kti(pt)L 2 \og(Mf 3 /M c )q sets the relevant mo- 
mentum scale. Thus R cb — > 1 more slowly for higher 
opacity. One can see this behavior reflected in the full 
numerical results shown in Fig. [2] for moderate suppres- 
sion, but that the extreme opacity q = 100 case deviates 
from Eq. Q. 

The maximum momentum for which string theoretic 
predictions for R cb can be trusted is not well understood. 
Eq. was derived assuming a constant heavy quark 
velocity. Supposing this is maintained by the presence 
of an electromagnetic field, the Born-Infeld action gives 
a "speed limit" of 7c = M 2 /A(T*) 2 [37]. The work of 
[19] relaxed the assumptions of infinite quark mass and 
constant velocity; nevertheless Eq. (JIJ well approximates 
the full results. Requiring a time- like endpoint on the 
probe brane for a constant velocity string representing a 
finite mass quark leads to [21] a parametrically similar 
cutoff, 



= 1 



2M 
VXT* 



AM 2 
\{T*f 



(5) 



There is no known limit yet for the dynamic velocity 
case. To get a sense of the pr scale where the AdS/CFT 
approximation may break down, we plot the momentum 
cutoffs from Eq. ^ for the given SYM input parameters 
corresponding to T*(to) and T*. These are depicted by 
"O" and "|" in the figures, respectively. 

Conclusions Possible strong coupling deviations from 
pQCD in nuclear collisions were studied based on a recent 
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AdS/CFT model of charm and bottom energy loss. The 
predicted nuclear modification factors, R-'aa' were found 
to decrease as a function of px, as compared to increasing 
as predicted from pQCD. The distinction between these 
dynamical models is amplified by studying the pr de- 
pendence of the double ratio R c aa (pt) / ' R\a(pt) > which 
clearly illustrates a quark mass independence of pQCD 
energy loss at asymptotically high momenta in contrast 
to the momentum independence and inverse mass de- 
pendence prediction of the AdS/CFT drag coefficient. A 
clear signal of novel nonpcrturbative physics would be 
the observation of a large deviation of the double ratio 
from unity at pr ^> Mb soon to be accessible at LHC and 
RHIC. 
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